Evaporative Cooling of Antiprotons to Cryogenic Temperatures 
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We report the application of evaporative cooling to clouds of trapped antiprotons, resulting in 
plasmas with measured temperature down to 9 K. We have modeled the evaporation process for 
charged particles using appropriate rate equations. Good agreement between experiment and theory 
is observed, permitting prediction of cooling efficiency in future experiments. The technique opens 
up new possibilities for cooling of trapped ions and is of particular interest in antiproton physics, 
where a precise CPT test on trapped antihydrogen is a long-standing goal. 



Historically, forced evaporative cooling has been suc- 
cessfully applied to trapped samples of neutral particles 
PQ, and remains the only route to achieve Bose- Einstein 
condensation in such systems [2j. However, the technique 
has only found limited applications for trapped ions (at 
temperatures ~ 100 eV [3]) and has never been realized 
in cold plasmas. Here we report the application of forced 
evaporative cooling to a dense (~ 10 6 cm~ 3 ) cloud of 
trapped antiprotons, resulting in temperatures as low as 
9 K; two orders of magnitude lower than any previously 
reported @]. 

The process of evaporation is driven by elastic colli- 
sions that scatter high energy particles out of the con- 
fining potential, thus decreasing the temperature of the 
remaining particles. For charged particles the process 
benefits from the long range nature of the Coulomb in- 
teraction, and compared to neutrals of similar density 
and temperature, the elastic collision rate is much higher, 
making cooling of much lower numbers and densities of 
particles feasible. In addition, intraspecies loss channels 
from inelastic collisions are non-existent. Strong coupling 
to the trapping fields makes precise control of the con- 
fining potential more critical for charged particles. Also, 
for plasmas, the self-fields can both reduce the collision 



rate through screening and change the effective depth of 
the confining potential. 

The ALPHA apparatus, which is designed with the 
intention of creating and trapping antihydrogen [5J, is 
located at the Antiproton Decelerator (AD) at CERN 
[BJ. It consists of a Penning-Malmberg trap for charged 
particles with an octupole-based magnetostatic trap for 
neutral atoms superimposed on the central region. For 
the work presented here, the magnetostatic trap was not 
energized and the evaporative cooling was performed in 
a homogeneous 1 T solenoidal field. 

FigureflJ, shows a schematic diagram of the apparatus, 
with only a subset of the 20.05 mm long and 22.275 mm 
radius, hollow cylindrical electrodes shown. The vacuum 
wall is cooled using liquid helium, and the measured elec- 
trode temperature is about 7 K. The magnetic field, indi- 
cated by the arrow, is directed along the axis of cylindri- 
cal symmetry and confines the antiprotons radially: due 
to conservation of angular momentum, antiprotons do 
not readily escape in directions transverse to the mag- 
netic field lines [7 . Parallel to the magnetic field, an- 
tiprotons are confined by electric fields generated by the 
electrodes. 

Also shown are the two diagnostic devices used in the 
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FIG. 1. a) Simplified schematic of the Penning-Malmberg 
trap used to confine the antiprotons and of the two diagnostic 
devices used. The direction of the magnetic field is indicated 
by the arrow, b) Potential wells used to confine the antipro- 
tons during the evaporative cooling ramp. The antiprotons 
are indicated at the bottom of the potential well (red), and 
the different wells are labelled by their on-axis depth. 

evaporative cooling experiments. To the left, antipro- 
tons can be released towards an aluminium foil, on which 
they annihilate. The annihilation products are detected 
by a set of plastic scintillators read out by photomul- 
tiplier tubes (PMT) with an efficiency of (25 ± 10)% 
per event. The background signal from cosmic rays is 
measured during each experimental cycle and is approx- 
imately 40 Hz. To the right, antiprotons can be released 
onto a microchannel plate (MCP) /phosphor screen as- 
sembly, allowing measurements of the antiproton cloud's 
spatial density profile, integrated along the axis of cylin- 
drical symmetry [9] . 

Before each evaporative cooling experiment a cloud of 
45,000 antiprotons with a radius of 0.6 mm and a den- 
sity of 7.5 x 10 6 cm~ 3 is prepared. The antiprotons are 
produced and slowed to 5.3 MeV in the AD, and as they 
enter our apparatus through the aluminium foil (see fig- 
ure [iji) they are further slowed. Inside the apparatus 
we capture 70,000 of the antiprotons in a 3 T magnetic 
field between two high voltage electrodes (not shown) 
excited to 4 kV [10]. Typically 65% of the captured an- 
tiprotons spatially overlap a 0.5 mm radius, pre-loaded, 
electron plasma with 15 million particles. The electrons 
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FIG. 2. The number of antiprotons lost from the well as it is 
lowered is integrated and plotted versus the well depth; note 
that, in time, the potentials are ramped from right to left in 
the figure. The measured number is corrected for the 25% de- 
tection efficiency. The curves are labelled in decreasing order 
of the temperatures extracted from an exponential fit, shown 
as the solid lines. The temperatures (corrected as described 
in the text) are: A: 1,040 K B: 325 K C: 57 K D: 23 K E: 19 K 
F: 9 K. As the antiprotons get colder, fewer can be used to 
determine their temperature, an effect described in Ref. [13] . 

are self-cooled by cyclotron radiation and in turn cool the 
antiprotons through collisions 0J. Antiprotons which do 
not overlap the electron plasma remain energetic and are 
lost when the high voltage is lowered. 

The combined antiproton and electron plasma is then 
compressed radially using an azimuthally segmented elec- 
trode (not shown) to apply a "rotating- wall" electric 
field, thus increasing the antiproton density [0J EJ H2] • 
The magnetic field is then reduced to 1 T and the parti- 
cles moved to a region where a set of low-noise amplifiers 
is used to drive the confining electrodes. Pulsed electric 
fields are employed to selectively remove the electrons, 
which, being the lighter species, escape more easily. The 
antiprotons remain in a potential well of depth 1500 mV 
(see Figure [Tja) . The quoted depth is the on-axis value 
due to the confining electrodes only. Space charge poten- 
tials are considered separately. This convention is used 
throughout the article. 

The temperature of the antiprotons can be determined 
by a destructive measurement in which one side of the 
confining potential is lowered and the antiprotons re- 
leased (in a few ms) onto the aluminium foil. Assum- 
ing that the particle cloud is in thermal equilibrium, the 
particles that are initially released originate from the ex- 
ponential tail of a Boltzmann distribution |13j , so that a 
fit can be used to determine the temperature of the parti- 
cles. Figure [2] shows six examples of measured antiproton 
energy distributions. 

The raw temperature fits in Figure [2] are corrected by 
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a factor determined by particle-in-cell (PIC) simulations 
of the antiprotons being released from the confining po- 
tential. The simulations include the effect of the time 
dependent vacuum potentials and plasma self- fields, the 
possibility of evaporation, and energy exchange between 
the different translational degrees of freedom. The sim- 
ulations suggest that the temperature determined from 
the fit is ~16% higher than the true temperature. Note 
that the PIC-based correction has been applied to all 
temperatures reported in this paper. The distribution 
labelled A in Figure [2] yields a corrected temperature of 
(1040 ± 45) K before evaporative cooling; the others are 
examples of evaporatively cooled antiprotons achieved as 
described below. 

To perform evaporative cooling, the depth of the ini- 
tially 1500 mV deep confining well was reduced by lin- 
early ramping the voltage applied to one of the electrodes 
to one of six different predetermined values (see examples 
on Figure [TJd). Then the antiprotons were allowed to re- 
equilibrate for 10 s before being ejected to measure their 
temperature and remaining number. The shallowest well 
investigated had a depth of (10 ± 4) mV. Since only one 
side of the confining potential is lowered, the evaporat- 
ing antiprotons are guided by the magnetic field onto 
the aluminium foil, where they annihilate. Monitoring 
the annihilation signal allows us to calculate the num- 
ber of antiprotons remaining at any time by summing all 
antiproton losses and subtracting the measured cosmic 
background. 

Figure[3^, shows the temperature obtained during evap- 
orative cooling as a function of the well depth. We ob- 
serve an almost linear relationship, and in the case of 
the most shallow well, we estimate the temperature to 
be (9 ± 4) K. The fraction of antiprotons remaining at 
the various well depths is shown on Figure [3]d, where it 
is found that (6 ± 1)% of the initial 45,000 antiprotons 
remain in the shallowest well. 

We investigated various times (300 s, 100 s, 30 s, 
10 s, 1 s) for ramping down the confining potential from 
1500 mV to 10 mV. The final temperature and fraction 
remaining were essentially independent of this time ex- 
cept for the 1 s case, for which only 0.1% of the particles 
survived. 

A second set of measurements was carried out to deter- 
mine the transverse antiproton density profile as a func- 
tion of well depth. For these studies the antiprotons were 
released onto the MCP/phosphor screen assembly (see 
Figure [ij,), and the measured line-integrated density pro- 
file was used to solve the Poisson-Boltzmann equations 
to obtain the full three dimensional density distribution 
and electric potential [14] , 

A striking feature of the antiproton images was the ra- 
dial expansion of the cloud with decreasing well depth, 
from an initial radius r$ of 0.6 mm to approximately 
3 mm for the shallowest well. If one assumes that all 
evaporating antiprotons are lost on the axis, where the 
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FIG. 3. a) Temperature vs. the on-axis well depth. The 
error is the combined statistical uncertainty from the tem- 
perature fit and an uncertainty associated with the applied 
potentials (one a). The model calculation described in the 
text is shown as a line, b) The fraction of antiprotons re- 
maining after evaporative cooling vs. on-axis well depth. The 
uncertainty on each point is propagated from the counting 
error (one a). The initial number of antiprotons was approx- 
imately 45,000, for an on-axis well depth of (1484 ± 14) mV. 



confining electric field is weakest, no angular momentum 
is carried away in the loss process. Conservation of total 
canonical angular momentum [7] would then predict that 
the radial expansion of the density profile will follow the 
expression N /N — (r 2 ) / (tq) when angular momentum 
is redistributed among fewer particles into a new equilib- 
rium. Here Nq is the initial number of antiprotons and N 
and r are, respectively the number and radius after evap- 
orative cooling. We find that this simple model describes 
the data reasonably well. 

To predict the effect of evaporative cooling in our trap 
we modelled the process by solving the rate equations 
describing the time evolution of the temperature T, and 
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FIG. 4. The evaporation timescale vs. applied on-axis well 
depth for five experiments with different ramp times. The 
on-axis well depth at the end of the ramp is (10 ± 4) mV 
in all cases. 7 -1 is indicated by the dot-dashed line, model 
calculations are shown as lines and, dashed lines indicate 77 < 
4 in the model calculation. 



the number of particles trapped N [TS] : 

dN _ N 
dt r e . 

dT T 

— — = —a 



P. 



(1) 
(2) 



where T ev is the evaporation timescale, and a the excess 
energy removed per evaporating particle. A loss term 7 
(1 x s _1 per antiproton) was added to allow for the 
measured rate of antiproton annihilation on the residual 
gas in the trap, and a heating term P was also included 
to prevent the predicted temperatures from falling below 
the measured limits. The value of P is determined by 
calculating Joule heating from the release of potential 
energy during the observed radial expansion. Assuming 
that the expansion is only due to particle loss and the 
conservation of angular momentum, as described above, 
we find P to be of order (—dN/dt)x5 mK. 

Following the notation in Ref. [TS] we calculate a as: 



77 + k 



- 1, 



5 + 3/2 " (3) 

where 77 is the potential barrier height, k the excess ki- 
netic energy of an evaporating particle and, 5 + 3/2 the 
average potential and kinetic energy of the confined par- 
ticles. The variables 77, k and 6 are the respective ener- 
gies divided by fcgT, where ks is Boltzmann's constant. 
Ignoring the plasma self-potential we let 8 = 1/2 and, 
following [T5] set k = 1. 

Using the principle of detailed balance, r ev can be re- 
lated to the relaxation time r co / for antiproton-antiproton 
collisions with: 

Tev _ \/2 
Tcol 3 



-77 e ' 



(4) 



being the appropriate expression for one dimensional 
evaporation (16) . Note that this expression is an approx- 
imation only valid for 77 > 4 [16] and, that the ratio is a 
factor of 77 higher than the corresponding expression for 
neutrals [T5]. For t co i we use the expression calculated 
in Ref. [T7]. 

From the antiproton density profiles, we estimate the 
central self-potential to be about 1.5 fiY per antiproton. 
In the initial 1500 mV well, this can be ignored, but with 
e.g. 20,000 antiprotons left in a 108 mV well a total self- 
potential of 30 mV changes our estimate of rj from about 
19 to 14 by reducing the energy required to evaporate 
from the well. The corresponding reduction in r ev is two 
orders of magnitude. Thus, the effect was included in the 
model, decreasing the predicted value of N by almost an 
order of magnitude in some cases. 

Figure [4] shows a comparison between the measured 
and the predicted evaporation timescales, calculated as 
r ev — — {d\ogN/dt + 7) _1 , for five different ramp times 
to the most shallow well. The good agreement between 
the measurement and the model calculations shows that 
we can predict r ev over a wide range of well depth and 
ramp times, giving good estimates of N in equation ([I]). 
The magnitude of 7 _1 is indicated on the figure, showing 
that antiproton loss to annihilations on residual gas is a 
small effect. 

In the calculations, shown as lines on Figure [3j we 
have modelled the temperature and remaining number 
as a function of the well depth. Of most interest is panel 
a), where reasonable agreement between the measured 
and the predicted temperatures validates the choice of 
a and P. With the measurements available we cannot 
exclude, but only limit, other contributions to P to be of 
similar magnitude as that from expansion-driven Joule 
heating. One such contribution could be electronic noise 
from the electrode driving circuits coupling to the an- 
tiprotons. The model also explains the excessive particle 
loss observed in the 1 s ramp. Here r ev becomes too short 
for rethermalization and 77 is forced below one. 

Evaporative cooling represents a strong candidate to 
replace electron cooling of antiprotons as the final cool- 
ing step when preparing pure low temperature samples 
of antiprotons for production of cold antihydrogen. Such 
samples can be used either as a target to create antihy- 
drogen through charge exchange with positronium atoms 
[181 119] or having no intrinsic cooling mechanism they 
can be precisely manipulated to have a well defined en- 
ergy relative to an adjacent positron plasma into which 
they are later injected. The technique eliminates electron 
related difficulties, such as heating by plasma instabili- 
ties [20] and centrifugal separation 21] in mixed plas- 
mas, absorption of thermal radiation and high frequency 
noise from the electrode driving circuit, and pulsed-field 
heating during removal. In principle, temperatures lower 
than the temperature of the surrounding apparatus are 
obtainable. Perhaps even sub-Kelvin plasmas, foreseen 
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to be a prerequisite for proposed measurements of grav- 
itational forces on antimatter [22j . can be achieved. 

Despite reducing the total antiproton number, we have 
increased the absolute number of antiprotons with an en- 
ergy below our 0.5 K neutral atom trap depth by about 
two orders of magnitude compared to the initial distri- 
bution; from less than one to more than 10 in each ex- 
periment. This greatly improves the probability of pro- 
ducing trappable antihydrogen in many proposed antihy- 
drogen production schemes [TH1 HH1 1231 HI] ■ Preliminary 
studies indicate that repeating the experiment with our 
octupole-based atom trap energized does not change the 
outcome. Electronic noise in the electrode driving circuit 
and the ability to precisely set the confining potential will 
most likely limit the lowest temperature achievable in any 
given apparatus. 

The evaporative cooling of a cloud of antiprotons down 
to temperatures below 10 K has been demonstrated. In 
general, the technique could be used to create cold, non- 
neutral plasmas of particle species that cannot be laser 
cooled. 
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